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Abstract 



CJ , The response of laser interferometers to gravitational waves has been calculated 

in a number of different ways, particularly in the transverse-traceless and the local 

Vh . Lorentz gauges. At first sight, it would appear that these calculations lead to 

different results when the separation between the test masses becomes comparable to 
the wavelength of the gravitational wave. In this paper this discrepancy is resolved. 
We describe the response of free test masses to plane gravitational waves in the 
coordinate frame of a local observer and show that it acquires contributions from 
three different effects: the displacement of the test masses, the apparent change 
in the photon velocity, and the variation in the clock speed of the local observer, 
all of which are induced by the gravitational wave. Only when taken together do 
these three effects represent a quantity which is translationally invariant. This 
translationally-invariant quantity is identical to the response function calculated in 
the transverse-traceless gauge. We thus resolve the well-known discrepancy between 
the two coordinates systems, and show that the results found in the coordinate frame 
of a local observer are valid for large separation between the masses. 
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1 Introduction 

Searches for gravitational waves are now conducted with laser interferometers in which 
the test masses for sensing gravitational waves are separated by distances of several 
kilometers PJE]- Variations in the proper distance between these test masses, which may 
be caused by gravitational waves, are measured with light. The response of the laser 
interferometers to gravitational waves is usually calculated in the transverse-traceless 
(TT) gauge [3j- The main assumption for all such calculations is that the test masses 
of the laser interferometers are inertial, i.e. accelerations of the test masses in the 
direction of the probe laser beam are negligible. A substantial engineering effort has 
been made to meet this requirement. Placed in ultra-high vacuum and isolated from 
the ground by multi-layer stacks and actively-controlled suspensions, these test masses 
become practically inertial at frequencies far above the suspension resonances. 

A significant change in the attitude toward the laser interferometers took place with 
the introduction of optical springs in the last few years [1J03III]. The optical spring is 
produced by the pressure of light on the test masses, which in the advanced interfer- 
ometer configurations can lead to amplification of the gravitational wave signal. For 
this amplification the resonance frequency of the optical spring must be matched with 
the frequency of the expected gravitational waves. In this case, the main assumption of 
the TT-gauge - the requirement of test mass inertiality - can no longer be made. This 
problem can be overcome if one uses the coordinates of a local observer for which no 
requirement of test mass inertiality is needed. Although this coordinate system has been 
frequently used to describe the response of resonant bar detectors jZj, its application to 
laser interferometers thus far has been only occasional. 

The coordinates of a local observer, also known as the local Lorentz gauge, have 
a long history. Comparison of the TT coordinates and the coordinates of a local ob- 
server is given in [2j with the curious observation that they yield different answers for 
geodesic deviation when separation between the geodesies becomes comparable to the 
wavelength of the gravitational wave (Exercise 37.6 in [2]). As a result, the coordinates of 
a local observer were considered not suitable for large separation between the geodesies. 
Nonetheless, the studies of the effects of gravitational waves in the coordinates of a local 
observer continued and over the years led to a number of interesting results. Several 
insightful papers have been written about the role of the coordinate system in the detec- 
tion of gravitational waves |HliniE3^B- Some of the calculations in these papers rely on 
the Fermi normal expansion as a means to build the coordinate frame of a local observer. 
Explicit transformations from the TT-coordinates to the coordinates of a local observer 
have been constructed and analyzed |121 1131 HI] . The response of the interferometric 
gravitational wave detector calculated in the TT-coordinates was transformed into the 
coordinates of the local observer H^J Hj)| . 

Despite of all these efforts, the coordinates of a local observer have remained an 
obscure gauge even to this day. One of the reasons for the lack of understanding is the 
avoidance of the local Lorentz gauge which is largely influenced by the disagreement 



between this gauge and the TT-gauge. Another drawback associated with the coordi- 
nates of a local observer is the lack of consistent mathematical formalism. To derive any 
nontrivial result in these coordinates, one usually starts with TT-gauge and then obtains 
the answer by complicated coordinate transformations. In this paper we show that it 
is possible to calculate the effects of the gravitational wave directly in the coordinates 
of a local observer. There is no need to start with the TT-gauge and no need to use 
the transformation rules to go from one coordinate frame to the other. In particular, 
we describe the response of test masses to gravitational waves when their separation 
is comparable to or greater than the wavelength of the gravitational wave. In this ap- 
proach, several effects must be combined to obtain a consistent test mass response. In 
the end, the discrepancy between the two gauges is resolved. 

The presentation in this paper is such that only a few concepts from differential 
geometry are used. Often, abstract mathematical derivations are replaced with those 
based on simple physical arguments, and many formulas are deliberately presented in 
the Newtonian form after they have been derived in general relativity. The motivation 
for this approach is two-fold. On one hand, it allows us to focus on physics of the 
problem and set aside mathematical details which can be overwhelming. On the other 
hand, such an approach allows us to assume the standpoint of a "Newtonian physicist" 
PI) conducting experiments in a laboratory environment and describing the outcomes 
of these experiments in the familiar Newtonian terms, even though they represent the 
effects in general relativity. 

2 The Coordinates of Transverse Traceless Gauge 

We begin with a brief overview of the TT-gauge. This digression will allow us to in- 
troduce the test mass response function which will be needed later for comparison. 
Subsequent calculations, however, do not rely on the TT-gauge in any way. 

In the TT-gauge the metric which describes a plane polarized gravitational wave 
propagating in flat space-time is given by 



QlLV 



/-I \ 

1+h 

1-h 



(1) 



where h = h(t + z/c) represents the amplitude of the "+" polarization 3. For all 
anticipated astrophysical sources, the amplitude of gravitational waves upon their arrival 
to Earth is expected to be extremely small, typically \h\ ~ 10~ 21 or less, which justifies 
the use of the perturbation method in the following calculations. 

A special property of the TT-coordinates is that an inertial test mass, which is 
initially at rest in these coordinates, remains at rest throughout the entire passage of 
the gravitational wave JS1E1- Here, the use of words "at rest" requires clarification: 
they only mean that the coordinates of the test mass do not change in the presence of 




Figure 1: Bouncing photon in the coordinates of the TT-gauge. 



the gravitational wave. The proper distance between any two test masses changes even 
though their coordinates remain the same. A convenient way to analyze variations in 
the proper distance is by means of "bouncing photons" ^7j- For example, a photon can 
be launched from one test mass to be bounced back by the other, as shown in Fig.^ For 
simplicity we assume that the test masses are located along the x-axis of the coordinate 
system. In this case, the interval takes the form: 



ds 2 = -c 2 dt 2 + [1 + h(t)] dx 2 



(2) 



The condition for a null trajectory (ds = 0) gives us the coordinate velocity of the 
photon: 



j_fdx\ 2 _ c 2 



dt) ~ 1 + h(t) ' (3) 

which is a convenient quantity for calculations of the photon propagation times between 
the test masses. As we know, the coordinates of the test masses, x a = I and Xf, = I + L, 
do not change under the influence of gravitational wave. Therefore, the duration of the 
forward trip can be found as 
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dx 
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where t' = t — (I + L — x)/c. To first order in h, this integral can be approximated as 

l+L 



T 1 (t)=T + ±- J h(t')dx, 



(5) 



where T = L/c is the light transit time in the absence of gravitational waves. Similarly, 
the duration of the return trip would be 

l 
T 2 (t) = r + 1 J h(t')(-dx), (6) 

l+L 

though now the retardation time is given by t' = t — (x — l)/c. 

The round trip time can then be found by adding T2(t) and T\[t — T2(i)]. The latter 
can be approximated by T\(t — T) because the difference between the exact and the 
approximate values is second order in h. Therefore, to first order in h, the duration of 
the round trip can be defined as 

T T . t .(t)=T l (t-T) + T 2 (t). (7) 

Deviations of this round-trip time from its unperturbed value, 2T, are then given by 

l+L _ 

dx. (8) 
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Even though I explicitly enters this equation, 5T does not depend on I. This observation 
implies that the choice of the origin for this coordinate system does not affect 5T; in 
other words, the result is translationally invariant. 

The deviation in the round-trip time, Eq.(JHJ), can also be written in the Laplace or 
Fourier domain. Laplace transformations are commonly used to analyze linear responses 
of interferometric gravitational- wave detectors ^H], and sometimes are easier to inter- 
pret than their time-domain equivalents. Define the Laplace transform of an arbitrary 
function of time h{t) by 

oo 

h(s) = f e~ st h(t) dt. (9) 

o 
Then the Laplace-domain version of Eq. © can be written as 

S JM = C( S )h(s), (10) 

where C(s) represents the response of test masses to gravitational waves: 

1 - e~ 2sT 

C(s) = - — . (11) 

w 2sT y J 

A number of derivations of this result, some quite different from ours, can be found in 
literature, for example in [191 I2U1 I2T1 \2'2\ and more recently in |231I18| . 

There are several reasons why the above picture is not satisfactory from a physical 
point of view, even though it is mathematically sound. The main problem with the 



coordinates of the TT- gauge is that they can hardly be realized in the experiment. In 
fact, they cannot be implemented in the laboratory environment on Earth because the 
coordinate grid of the TT-frame must be changing in unison with the passing gravi- 
tational wave, the effect commonly known as "breathing of the frame." (They may, 
however, be realized in space with a network of freely-falling satellites.) Consequently, 
the application of the above calculations to ground-based gravitational-wave detectors 
becomes problematic. For physicists working with these detectors, it is sometimes not 
clear how the results derived in the TT-coordinates can be used in experiments when 
these coordinates are not available in practice. 

Another problem, which is closely related to the previous one, comes from the as- 
sumption of test mass inertiality. Namely, the above derivation of the photon round-trip 
time was based on the premise that the coordinates of the test masses do not change 
under the influence of the gravitational wave, an assumption which is true only when the 
test masses are inertial. The test masses in laser gravitational-wave detectors constantly 
undergo accelerations in response to various forces acting on them and thus are never 
truly inertial. One can argue that these accelerations typically occur at frequencies of 
the suspension resonances which are well below the frequencies of the anticipated gravi- 
tational waves. However, in advanced interferometer configurations the accelerations of 
test masses will also be caused by the radiation-pressure variations which are intended 
to occur at the frequency of anticipated gravitational waves. Therefore, the assumption 
of test mass inertiality, which is most effective in the TT-gauge, becomes too restrictive 
for more realistic calculations. These problems do not occur if one uses the coordinates 
of a local observer. 

3 The Coordinates of Local Observer 

An observer in a laboratory environment on Earth typically uses the coordinate system 
in which the space-time is locally flat ^7j, and the distance between any two points 
is given simply by the difference in their coordinates in the usual sense of Newtonian 
physics j3j. In this reference frame, gravitational waves manifest themselves through the 
tidal forces which they exert on the masses j2j- To describe the tidal forces we consider a 
test mass which is free to move in the horizontal plane (z = 0). For simplicity, we assume 
that this plane coincides with the wavefront of the gravitational wave, and that the x 
and y directions of the coordinate system match the polarization of the gravitational 
wave. Then the tidal acceleration of the test mass caused by the gravitational wave [H] 
is given by 

1 ■■ 
x = +-hx, (12) 

v = -\hy- (13) 



Equivalently [HJE1[2I], one can sa Y that there is a gravitational potential: 

^r,t) = -\h(t)(x 2 -y 2 ), (14) 

which generates the tidal forces, and that the motion of the test mass is governed by 
the Newton equation: 

r = -V$. (15) 

The potential is not static, and therefore the energy of the test mass is not conserved. 
This is the Newtonian version of the theorem from general relativity which states that 
gravitational waves must supply energy to test masses to become detectable in experi- 
ments |25j . 

In the post-Newtonian approach to general relativity (see also Appendix 0), the 
gravitational potential is related to the time component of the metric: 

500 = -1 - \ *. (16) 

In what follows we will frequently use a perturbation expansion, keeping only terms 
which are first order in h, and therefore rely on the assumption that |^|/c 2 <1. To 
satisfy this condition, we require that the spatial coordinates x and y do not extend 
indefinitely. This limitation, however, will not restrict us in any way. Indeed, for 
gravitational waves with the largest expected amplitudes (\h\ ~ 1CP 21 ) and the highest 
detectable frequencies (~ 10 kHz), the restriction on the spatial coordinates implies that 
|x|, \y\ <C 10 14 m, which is always satisfied in the laboratory environment on Earth. 

The solution to Eqs. lfT2|) - lji;-{JI is usually found using the perturbation method j3j. 
To first order in h, the displacements of the test mass caused by the gravitational waves 
are given by 

5x(t) = +^x h(t), (17) 

Sy(t) = ~yoh{t), (18) 

where xq and yo are the initial (unperturbed) coordinates of the test mass. This notion 
is regarded as the major difference between the coordinates of a local observer and 
the coordinates of the TT-gauge, in which the test masses were not moving under the 
influence of the gravitational wave. 

4 Requirement of Translational Invar iance 

An interesting feature of the local Lorentz gauge is the coordinate dependence of the 
tidal forces - they can be changed by a mere shift of the origin of the coordinate system: 

x — > x + X, and y — ► y + Y. (19) 



The same applies to the test mass displacements, Eas.(|17|l and (|18j). This is the earliest 
indication that the coordinates of a local observer are not as simple as they may seem. 
However, at this point, the coordinate dependence seems to be quite harmless, and 
we can entertain the notion that it can be removed simply by considering the relative 
motion of test masses. 

As before, we probe the geometry of space-time with a bouncing photon. Consider 
two test masses with coordinates x a and Xb and assume that the photon is launched 
from one test mass and is bounced by the other. Let the unperturbed values for the test 
mass coordinates be 

x a = I, and Xb = I + L, (20) 

and the unperturbed propagation time between the masses be 

T=-. (21) 

c 

From Ea. ()17|) we find that the displacements of the test masses under the influence of 
the gravitational wave are 

5x a (t) = ~lh(t), (22) 

Sx b (t) = ^(l + L)h(t). (23) 

Then the relative displacement, commonly defined as 

6L(t) = Sx b (t) - 5x a (t) 

= \Lh(t), (24) 

would obviously be independent of / and therefore independent of the choice of the origin 
for these coordinates, as we anticipated. Equation PH) . often written as 

^-± am. (*> 

is widely used to describe the strain induced by gravitational waves on bar detec- 
tors. However, its application to laser interferometers immediately runs into a problem. 
Namely, the change in the round-trip time calculated from Eq.(|25|) would be 

^ = SM, (26) 

which is different from the one obtained in the TT-gauge, Eq. (|lU|) . This is the precise 
origin of the well-known discrepancy between the two coordinate systems. One of the 
earliest accounts of this discrepancy appears in Exercise 37.6 of [3], which also suggests 
that the correct answer for the photon propagation time must be obtained in the coordi- 
nates of the TT-gauge. It is sometimes assumed that the discrepancy occurred because 



of the application of Eas. (|12|) - ()13|) beyond their limits of validity. The actual cause of 
the discrepancy lies in the neglect of the effects of gravitational redshift, as will be shown 
below. 

Historically, the discrepancy was not viewed as a serious problem when the searches 
for gravitational waves were conducted with bar detectors. The relatively small size of 
a bar detector (a few meters) implies small separation for its constituent parts, in which 
case the difference between the two coordinate systems becomes negligible. Indeed, for 
gravitational waves with wavelength much greater than the separation between the test 
masses, \sT\ <C 1 and therefore C(s) rj 1, which makes Eq. ljlOJI equivalent to Eq. (|26j) . 
The situation became rather different with the arrival of long-baseline laser interferom- 
eters. In these detectors the test masses for sensing gravitational waves are separated 
by distances of several kilometers, and the long- wavelength approximation, \sT\ <C 1, 
becomes hard to justify. Furthermore, recent studies [Sfj] have shown that these inter- 
ferometers are capable of detecting gravitational waves with wavelengths comparable to 
their arm-length, \sT\ ~ 1, thus operating entirely outside the long-wavelength regime. 

5 Requirement of Causality 

For large separation between the test masses, the definition for relative displacement, 
Eq. (|24[) . becomes unphysical. In this definition the two test masses are taken at the 
same time and therefore cannot be in causal connection. The definitions for the relative 
test-mass displacement which are appropriate for the bouncing photon can be written 

as 

<5Li(£) = fa 6 (*)-&c (f-Ti), (27) 

SL 2 (t) = 5x b (t-T 2 )-5x a (t), (28) 

where T\ and T 2 are the photon propagation times for the forward and return trip 
correspondingly. According to these definitions, the displacement of one test mass is 
compared with the displacement of the other at a later time to allow for finite delay 
from the light propagation, as can be seen from Fig. [21 Note that the propagation times 
T\ and T 2 in Eas. (|2"7j) - (|28|) can be replaced with their nominal value T because the test 
mass displacements are already first order in h. 

The total change in the distance between the masses in one round-trip of light would 
be 

5L T . t .(t) = SL^t-^ + SL^t) 

= 2 5x b (t-T)-6x a (t)-6x a (t-2T). (29) 

An explicit formula for this length change written in terms of the amplitude of the 
gravitational wave is 

5L T . t .(t) = (I + L)h(t - T) - - I h(t) - - I h(t - 2T). (30) 




Figure 2: Bouncing photon in the coordinates of a local observer. 



This quantity represents causal length variations in laser interferometers for gravitational 
wave detection. Note that <5L r .t. is not translationally invariant, despite the fact that it 
represents the relative displacement of the test masses. This is the price one has to pay 
for satisfying the causality condition. 

Changes in the distance, Eq. pUj) . lead to changes in the round-trip time for photons 
propagating between the masses: 



S x T(t) 



h(t - T) - n [h(t) - 2h(t -T) + h(t - 2T)\ , 



where we introduced a dimensionless parameter 

I 



" = 2L' 



(31) 



(32) 



The presence of this parameter in the subsequent formulas will indicate the loss of 
translational invariance. The Laplace-domain version of Eo. pij) can be written in a 
manner similar to Eq. ()lU|) . namely 



5 x T(s) 



D x (s)h{s). 



where D x (s) is the corresponding response function 

D x (s) = e~ sT - ,*(! - e^ T ) : 



(33) 



(34) 



Note that D x (s) depends on the choice of the origin for this coordinate system. At 
first it may seem that this loss of translational invariance is natural. After all, the 
potential explicitly depends on coordinates, which in classical mechanics usually means 



that the symmetry with respect to translations is lost. However, such a conclusion would 
contradict our physical intuition which maintains that all locations on the wavefront of 
the plane gravitational wave must be equivalent. This implies that physical quantities 
must be the same no matter where on this plane they are measured, even though the 
potential explicitly discriminates between different locations on the plane. We will see 
shortly that this is indeed the case and that translational invariance is restored, but 
only when another significant effect is added to the picture: the gravitational redshift 
of light propagating between the masses. 

6 Distributed Gravitational Redshift 

We have calculated variations in the photon round-trip time which come from the motion 
of the test masses induced by the gravitational wave. In this calculation, we implicitly 
assumed that the propagation of the photon between the test masses is uniform, as if it 
were moving in flat space-time. However, the presence of the tidal forces indicates that 
space-time is curved. As a result, the bouncing photon will experience a gravitational 
redshift. There will be two such effects in the following calculations. The first will be 
called the distributed gravitational redshift because it requires spacial separation, the 
second will be called the localized gravitational redshift because it occurs at a single 
point in space. 

The distributed gravitational redshift can be calculated as follows. Consider the 
interval for photons propagating along the a;- axis: 

ds = goo c dt + dx , (35) 

where goo is the time component of the metric, Eq. ljffij) . The condition for a null trajec- 
tory (ds = 0) gives us the coordinate velocity of the photons: 

2 - f ^} 2 = c 2 + 2<S>(t,x). (36) 

To first order in h, the velocity can be approximated by 



v ~ ± c 






(37) 



where + and — correspond to the forward and return trip, respectively. 

Knowing the coordinate velocity of the photons, we can define the propagation time 
for the photon traveling between the masses: 

Xb(t) X a (t) 

T l (t)= I -, and T 2 (t)= f ^, (38) 

J v J v 

Xa(t~Tl) x b (t-T 2 ) 
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in accordance with Fig. [3 We will not attempt to calculate these integrals directly. 
Such calculations would be complicated because the boundaries of these integrals are 
changing with time: 



x a {t) = l + 5x a (t), 
x b (t) = l + L + 5x b (t). 



(39) 
(40) 



Fortunately, we do not need to calculate the contributions of the boundary terms to these 
integrals. To first order in h, these contributions can be approximated by 5L\(t)/c and 
8L/2(t)/c (see Eas.([27|) and (HE}). Therefore, the combined effect of the varying bound- 
aries is given by S X T, previously found in Eq. (j31j) . Thus, we only calculate the times for 
photon propagation between the fixed boundaries: I and I + L. Such propagation times 
will be denoted here by ATi^ to be distinguished from T\^- 

In the forward trip, the propagation time between the fixed boundaries is 



l+L 



ATi(t) 



dx 



v(t',x) 

l+L 



T -— I $(t',x)dx, 



(41) 
(42) 



where t' is the retardation time which corresponds to the unperturbed photon trajectory: 
t' = t — (I + L — x)/c. Similarly, the propagation time between the fixed boundaries in 
the return trip is 



AT 2 (i)=T-l J <f>(t',x)(-dx), 

l+L 



(43) 



though now the retardation time is given by t' = t — (x — l)/c. The round-trip time 
for photons traveling between the fixed boundaries can be found by adding AT 2 (i) and 
ATi(t — T). Deviations of this round-trip time from its unperturbed value, 2T, are given 
by 

l+L 



5 v T{t) = -1 



x — I \ ( x — I 
$ [t-2T + ,x) +$ [t ,x 



dx. 



(44) 



After replacing the potential with its explicit form, Eq. ()14(l . we obtain the formula for 
5 V T in terms of the amplitude of the gravitational wave: 



l+L 



5 vT{t) = -^ 



m x — 1\ ■■ ( x — I 
h[t-2T + )+h[t 



x 2 dx. 



(45) 



This quantity represents the effect of the distributed gravitational redshift 6. 
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Equation (|45|) bears close similarity with Eg. ([%]). as both formulas represent cumu- 
lative effects of the gravitational wave. However, unlike Eq.(jHJ), which is translationally 
invariant, Eq. ()45|) is not, as can be seen from the presence of x 2 -factor in the integrand. 
A better way to analyze the loss of translational invariance would be to rewrite the 
result in the Laplace domain: 

6 -^ = D v (s)h(s), (46) 

where D v (s) is the corresponding response function: 

o.i.) _ ± f ( } - e -^)- e -.T + 

fl ^-e-' T f + i i 2 (l-e-'" T )sT^ (47) 

The terms proportional to fi and /i 2 represent the dependence of the response function 
on the choice of the origin for this coordinate system. 

We can now combine the variations in the photon propagation time which are caused 
by the motion of the test masses with those caused by the distributed gravitational 
redshift. The resulting round-trip time would be 

T, t . = 2T + 5 X T + 8 V T. (48) 

To this point, the combined effect of the gravitational wave is given by the sum: 

D x {s) + D v (s) = (^ + !?&) (l - e- 2sT ) . (49) 

By adding the two response functions we cancel the terms proportional to /U. However, 
the term proportional to fi 2 remains. As will be shown next, this term is related to the 
localized gravitational redshift. 

7 Localized Gravitational Redshift 

The last contribution to the photon round-trip time is also related to the gravitational 
redshift, although it is somewhat different from the distributed effect described above. 
In the presence of gravitational waves the clocks at different places run differently. The 
rate (dt*) of the clock which is located at x is related to the rate (dt) of the clock at the 
origin by 

dt* 2 = -g 00 (t,x)dt 2 , (50) 

which is the proper time at this location. In the above calculation of the photon round- 
trip time, Eq. (|48j) . we implicitly assumed that the time is measured with the clock at 
the origin: x = 0. The photon trajectories, however, begin and end at the location of 
the first test mass, a finite distance (I) away from the origin. As a result, the readings of 

12 



time become dependent on this distance. To avoid this problem, we shall measure time 
with the clock located at x = I. For this clock, the round-trip time is different from T rt ., 
Ea. (|48|) . The presence of the time-dependent gravitational potential affects the reading 
of this clock, causing it to register the round-trip time as 

t 
T* t .(t) = J sf-g 00 (t>,l)dt' 

t-Tr.t. 

t 

~ Tr.t.(t) + ^ J $(t',l)dt'. (51) 

t-Tr.t. 

To first order in h, the variation of the round-trip time due to this effect can be estimated 

as 

t 
5 t T(t) « ^ / $(t',l)dt' (52) 

t-2T 

= -&[[h(t)-kt-zn]- (53) 

This contribution to the round-trip propagation time comes from the non-uniformity of 
time flow caused by the presence of the gravitational wave. It will be called here the 
localized gravitational redshift. In the Laplace domain it can be written as 

^ = AWAW, (54) 

where Dt(s) is the corresponding response function 

D t (s) = -fi 2 (l - e~ 2sT ) sT. (55) 



Addition of this response function to Eq. (j49j) will cancel the // 2 -terms, giving us a 
translationally invariant result. 

We can now conclude that the change in the round-trip time caused by the gravita- 
tional wave consists of three contributions: 

5T = 5 X T + 6 V T + S t T, (56) 

which come from displacement of the test masses, changes in the coordinate velocity 

of the photons and variations in the clock rate. The combined result of these effects is 

given by the sum: 

1 - e~ 2sT 
D x (s) + D v (s) + D t (s) = ^ T , (57) 
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which is translationally invariant. Furthermore, the sum gives us a response function 
which is identical to C(s), Eq. (|11|). which proves that the two coordinate systems indeed 
yield the same answer for the observable photon round-trip time. 

The requirement of translational invariance played a special role in the above anal- 
ysis. The coordinate transformations, Ea.H19[). are a particular case of transformations 
known as changes of the origin, which in general relativity are accomplished with the 
help of parallel transports [T2|. Following the Newtonian style of our presentation, we 
referred to these transformations as translations and assumed that they represent a 
symmetry. The origin of this symmetry is related to the planeness of the gravitational 
wave [25] . 

8 The Round-Trip Phase of Light 

In the above picture, we considered the bouncing photon as a particle, assuming that 
there is a beginning and an end for the photon round trips. In practice, measurements of 
photon propagation times are usually done with optical interferometry in which photons 
are represented by continuous electromagnetic waves. We shall therefore briefly describe 
how the above calculations can be modified to become applicable to continuous waves. 
To be specific, we assume that the light is represented by a plane monochromatic wave 
with frequency uj and wavenumber k. In the absence of gravitational waves, such a wave 
is given by exp[i(ujt^kx)}. Then the photon trajectory introduced above would describe 
advancement of a surface of constant phase, whereas the photon velocity becomes the 
phase velocity of the wave. In this approach, the quantity of interest would be the 
round-trip phase, or more precisely, its variation caused by the gravitational wave. 

The first contribution to the round-trip phase comes from the motion of the test 
masses: 

ip x = -k5L T . ti = -uj 5 X T, (58) 

where <5L r . t . represents variations in the distance between the test masses, Eq. ()3U|) . and 
S X T represents the corresponding time variations, Eq. (|Hl|) . The second contribution 
comes from the variations in the phase velocity of the wave: 

^k = ^[$dx. (59) 

c 

Here we give a brief derivation of this result based on simple physical arguments. (An- 
other derivation, based on the solution of the eikonal equation, is given in Appendix 

El) 

In the presence of gravitational waves, the frequency and wavenumber are no longer 
constant; they become functions of position and time: Q(x,t) and K(x,t). Then the 
dispersion relation for the electromagnetic wave would read 

Q 2 = v 2 K 2 , (60) 
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where v is the phase velocity of the light, defined in Ea. (|3l)|) . To first order in h, Ea.(|6U|l 

can be written as 

k 
n-cK = -$. (61) 

c 

For a plane electromagnetic wave moving in the positive x-direction, an infinitesimal 

phase change is given by (ildt — K dx). Then the accumulated phase change can be 

found by integrating this quantity along the trajectory of a given wavefront. In doing 

so, we would find that the accumulated phase change vanishes because dx/dt = 0,/K. 

This result is quite natural, as traveling with the wavefront means following the surface 

of constant phase for which no phase change ensues. However, we must remember that 

we are not interested in the absolute phase change along the photon trajectory. Rather, 

we are interested in the variation of this phase change with respect to the unperturbed 

wave. Such a phase variation would be given by 



t/j k = (ndt-Kdx), (62) 

c 

where C denotes the unperturbed trajectory: dx/dt = ±c. Taken along the unperturbed 
photon trajectory, such an integral would yield a non-zero answer, which is equivalent 
to Ea. (|59j) . Note that the integral over the photon trajectory, Eq, ()59|) . has already been 
calculated, see Eq. (|H|) . Therefore, the phase change can then be written as 

ip k = -u) 5 V T, (63) 

where 5 V T is the corresponding variation in the round-trip time. 

We now can add this phase change to the phase change produced by the motion of 
the test masses, Eq . (|58[1 . There is no need to worry about the difference between k and 
K for this part. The displacements of the test masses are first order in h, and therefore 
any correction to k would result in second order terms. Thus, the combined effect is 
given by 

^ = ^ x + ^ k = -uj(5 x T + 5 v T). (64) 

As we already know from Eas.([48|) and (|49|) . this quantity is not translationally invari- 
ant, which means that it cannot be observed in the experiment. The change in the phase 
shift V represents the difference of the round-trip phases for two different electromag- 
netic waves: with and without the gravitational wave. Such a phase change cannot be 
measured in the experiment as the two waves cannot exist in the same space-time. To 
form an observable quantity, we shall compare the phase change of the probe electro- 
magnetic wave with that of a reference wave. The natural reference is the source itself, 
and therefore we need to find the phase change of the source. 

In flat space-time, the phase of the source would simply be ujt, and the phase shift of 
the source 2uT. In the presence of gravitational waves, the phase of the source becomes 
cot* , where t* is the proper time at the location of the source. Then the phase shift of 
the source can be found as 

u) [t*(t) -t*(t- r r . t .)] =uT* t Xt). (65) 
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Therefore, the change in this phase shift which is produced by the gravitational wave is 

ip so = w S t T, (66) 

where StT is given by Eq.(|52|). 

We can now compare the phase change of the moving wave, Eq . (|64|) . with that of 
a static source, Eq. (|66j) . The difference between the phase change of the wavefront for 
the electromagnetic wave returning to the source and the phase change of the source at 
that moment is 

Sip = i> - ipso- (67) 

In the explicit form this phase difference is given by 

Sip = -u(S x T + 5 V T + 5 t T). (68) 

As we already know from Eas.(|56|) and (|57|). this phase is translationally invariant and 
therefore represents an observable quantity. It is not surprising that this phase deviation 
is related to the time deviation by a simple formula: 

Sip = -to ST. (69) 

This result could have been guessed from simple dimensional analysis and the require- 
ment of translational invariance. The above derivation serves to explain the physical 
meaning of the relative phase shift and its constituent parts. In short, the motion of 
the test masses and the distributed gravitational redshift appear now as the phase shift 
of the traveling wave, whereas the localized gravitational redshift appears as the phase 
shift of the static source. 

9 Concluding Remarks 

We have shown that the response of test masses to gravitational waves in the local 
Lorentz gauge acquires contributions from three different effects: the motion of the 
test masses and the distributed and localized gravitational redshifts. Only when taken 
together do these effects yield an observable quantity. The approach developed in this 
paper has allowed us to calculate physical quantities directly in the coordinates of the 
local observer. In these coordinates, the assumption of the test mass inertiality is not 
required, and various forces acting on the masses can be added at will. We have provided 
a consistent framework for doing calculations in the coordinate system which is more 
natural for ground-based laser gravitational-wave detectors than the TT-gauge. 

To simplify the calculations, we introduced the three effects of the gravitational wave 
in a step-by-step fashion. At each step, mathematical derivations took advantage of the 
previous step. In retrospect, it is clear that a more direct way of doing the calculations 
would be to start with an abstract definition: 

T* t .=Jdt*, (70) 
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and then to proceed with integration over the photon trajectory 

T* t . = J y/-goo(l,f)de = J y /-g QQ (l,t>) 1 £^. (71) 

c* 

In this approach the contour C* would represent the actual photon trajectory: dx/dt = 
±v(x,t) which extends to the actual test mass positions: I + 5x a and I + L + 5xb- By 
evaluating various terms in the contour integral to first order in h, one would reproduce 
the above three contributions to the round-trip time. Although this approach may seem 
different from the one described in this paper, the mathematical equations and their 
physical interpretations would be essentially the same. 
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A Coordinate and Metric Transformations 

For completeness, we present here the transformations from the TT-gauge to the gauge 

of a local observer. Denote the coordinates of a local observer by x^ 1 and the metric in 

these coordinates by g^ u . Also, denote the coordinates of the TT-gauge by x 11 and the 

corresponding metric by g a p. The components of the metric in the TT-gauge, Eq.Q, 

can be written as 

/ \ 

ft, 

-h 

\ / 
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9flV — V/MV T 



(72) 



where r]^ v = diag{— 1, 1, 1, 1} is the Minkowski metric, and h = hit + z/c). From 
general relativity we know that the coordinate transformations, x^ — > x^, induce the 
transformations of the metric: 



9iiv 



dx a dx$ 
dx^ dx h 



9af3- 



(73) 



By definition, g^ v must become the Minkowski metric at the origin, and all its derivatives 
must vanish at this point. There are a number of metrics which satisfy these conditions. 
Here we consider one such choice El El • It can be obtained with the coordinate 
transformations, which to first order in h, are given by 



t = t--^h(x 2 -y 2 ), 

x = x hx, 

2 

V = V+^hy, 

z = z + — h{x 2 -y 2 ). 
4c 



(74) 
(75) 
(76) 
(77) 



The corresponding metric tensor can be obtained by performing the induced transfor- 
mation, Ea. (J73|) . To first order in h, the result is 



9/iu Vfiu 
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(78) 



where $ is a function of the new coordinates 

1 
4 



<I> 



h(t + z/c) (x 2 - y 2 



(79) 



In the post-Newtonian approach, this function becomes the potential for the tidal forces, 
which for z = appeared in Eq.(|14|), 

It is interesting to note that although the transformation rules Eas. ()74j) - H77|) are 
approximate, the metric Eq.(|78|) is an exact solution of Einstein equations |27j . This 
metric is generally known as the plane-front solution for strong gravitational waves 
|28| I29j . Further discussion of the relationship between the metric of the local observer 
and the exact solution can be found in |24j . 

The propagation of an electromagnetic wave in curved space-time is described by 
the eikonal ^ 301, which satisfies the equation: 



9' 



fUJ 



dm dm 

dx^ dx v 



0, 



(80) 
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where g^ v stands for the contravariant metric tensor. Its components are given by 



^v 



„[tV 



^-$00 $ \ 



\ $ 0-$/ 



For light propagating along the x-axis, the eikonal equation becomes 

V~^ V [da) = \dx~) ' 

The solution of this equation is described next. 



(81) 



(82) 



B Solution to the Eikonal Equation 

The eikonal equation can reduced to a linear differential equation by taking the square 
root of both sides of Eq. (|82|) and by keeping only the terms which are first order in h: 



dot dx 



del) 



(83) 



where ± correspond to the wave propagation in the positive and negative x-directions. 
The large unperturbed value of the eikonal satisfies Eq. ()83|) in the absence of the grav- 
itational waves ($ = 0), and is given by tot =F kx up to an additive constant. Therefore, 
to first order in h, the solution of the eikonal equation can be found as 



vl/ 1 = uot - kx + kl + 5^i, 

* 2 = ujt + kx- k(l + 2L) +5^>2, 



(85 



where 5^1,2 are the first order perturbations. For convenience we introduce the light 
cone coordinates: 






(ct + aO/2, 

(ct - x)/2, 



(86 
(87 



in which the photon world- lines become collinear with the £ and r\ axes, as shown in 
Fig. |21 In these coordinates, the first order perturbations satisfy the equations: 



0_ 
Or] 



S^fK 



k 



$. 



These equations allow direct integration: 
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Figure 3: Propagation of the electromagnetic waves in the light cone coordinates. The 
source is located at x = I and the turning point at x = I + L. 



Mi(Z,v) 



*(?,v)d? + fi(v), 
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k 
<5* 2 (£,??) = ^J *(Z,T/)dr/ + f2(Z), 
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where fi(r/) and f 2 (0 are arbitrary at this point. Transforming back to the coordinates 
x and t, we obtain the solution: 



6Vi(x,t) = -^[$(x',t-- — -)<y + fi(x,t), 



(92) 



l+L / 

6V 2 (x,t) = -» / $hr',i + ^— - ) dx' + f 2 {x,t). 



(93) 



The function /i is defined by the value of the eikonal at the source for the electro- 
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magnetic wave: v &i(/,t) = ujt*. To first order in h, this value can be found as 



*i(M) 



uJ 



J yJ-goo(l,V)dt' 



(94) 
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l + -q>(l,t') 



dt' 



(95) 



which must be the same as ut + fi(l,t), according to the definition Ea. (|51|) . We thus 
find the function f\ at the location of the source. Knowing that f\ is a function of ct — re, 
we can extend the values of f\ from the source location to the entire xt-plane: 



h(x,t) = - j m,t')dt'. 





(96) 



The function f 2 is defined by the value of the eikonal ^i at the turning point. The 
continuity of the eikonal implies that 



6-$ l (l + L,t) = 5y 2 (l + L,t). 



(97) 



From this condition we can find f 2 at the turning point. Knowing that f 2 is a function 
of ct + x, we can extend the values of f 2 from the turning point to the entire xi-plane: 



t _ 2T+ ^zi 



f2(x,t) 



- j ^{l,t')dt' + 
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®(x',t-2T+ ) dx. 
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The phase shift acquired by the electromagnetic wave in one round trip is given 
by the difference between the value of the eikonal at the beginning and the end of the 
propagation. To first order in h, this phase shift is given by 



Mt) = ^2{l,t)-5^ 1 {l ) t-2T). 
Simple algebra shows that this definition leads to 

l+L 

rp k (t) = -^ I $(x,t-- ) dx + 

I 

l+L 

-s f <S>(x,t-2T+?—— ) d.r. 
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which is an extended form of Ea. (|59|) . 

Finally, we give explicit formulas for Q and K in terms of the gravitational potential. 
These two quantities can be derived from the eikonal: 



n 






and K = q= 



dx ' 



(101) 



where — corresponds to the forward trip and + to the return trip. For example, in the 
forward propagation 



n(x,t) = u + -r(x,t) + -$(l,t-- 

c c \ c 

k „, k ( x — I 

K(x,t) - k + —r(x,t) + — <5>(l,t 



-j*(a;,t), 



cr cr V c / c* 

where T represents the non-stationary effect of the gravitational redshift: 



T(x,t) 



X 

i r d 



cj dt 
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x — X 



dx . 



Note that K can also be written as 

1 k 

K(x,t) = -Sl(x,t)--~$(x,t), 
c <r 

which leads directly to the dispersion relation, Eq. ljHTj) . 
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